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I. INTRODUCTION 


A. BACKGROUND 

Since the destructive effect of near proximity underwater explosions from 
torpedoes and mines was first documented in World War I, the U.S. Navy has pursued an 
experimental program to study these events and their effects. The Bureau of Ships 
organized the first concerted effort at the Norfolk Naval Shipyard in the 1930's. The 
introduction of more sophisticated weapons during World War II intensified and 
broadened research in this area. Since these early days the Navy has supported an 
experimental program to investigate methods for improving the resistance of ships and 


submarines to the effects of UNDEX. [Ref. 1] 


This effort has led the Navy to develop guidelines and specifications for the shock 
testing and hardening of shipboard equipment and systems. NAVSEA 0908-LP-000- 
3010A [Ref. 2] and MIL-S-901D [Ref. 3] are examples of this guidance. The shock 
resistance validation is then conducted through shock trials as required in OPNAVINST 
9072.2 [Ref. 4]. 


OPNAVINST 9072.2 mandates a shock trial for the first ship of every shock 
hardened class of surface ships. Testing the lead ship enables LFT&E to fulfill its role in 
correcting design deficiencies early on. These tests, although conducted under conditions 
more benign than realistic threats due to cost and crew safety concerns, provide insight 
into surface ship vulnerabihties to close proximity underwater bursts. The test subjects 
the ship to a series of UNDEX shocks generated by explosive charges at large standoff 
distances. The resulting uniform side-on shock is far from representative of realistic 
threats encountered in today's tactical environment. Additionally, the test 1s setup to 


eliminate the "whipping" response of the ship's structure due to the oscillations of the 


UNDEX bubble. [Ref. 5] 


Successful planning and completion of ship shock trials can be time consuming 
and expensive. The June 1994 shock trial of the USS John Paul Jones (DDG 53) began 


four years prior to the actual test date, cost tens of millions of dollars to complete and was 


met with lawsuits filed against the Navy by environmentalist groups [Ref. 6]. These cost, 
safety and environmental concerns are leading to greater interest in computer modeling 


and simulation of ship shock trials. 


Advances in computer technology, finite element modeling and simulation have 
lead to the question: "Can computer simulations replace ship shock trails?" The current 
state of computer technology is inadequate to reliably simulate the complex situation 
under investigation in the area of UNDEX. This analysis involves not only modeling the 
dynamic response of a very large complex structure but also the fluid-structure 
interaction of the UNDEX shock wave and the ship's hull. Reliable data from 
simulations would require high fidelity, full ship structural models that include wiring 
and piping systems as well as structural components like beams and shell plating. All 
these items contribute to the damping of the ship structure. Some amount of the 
surrounding fluid must also be modeled to capture the effects of the cavitation caused by 
the UNDEX shock waves. Additionally, the size of the finite elements and the time-step 
utilized in the simulation define the minimum length and time scales of the responses that 
can be captured. Shock propagation through the structure occurs in microseconds, while 
a global structural response, like whipping, takes seconds. The enormous model size, and 
wide range of length and time scales, required to fully capture the UNDEX test would 


demand tremendous computational resources. [Ref. 5] 


Despite these limitations modeling and simulation of ship shock trials can play a 
significant role in the ship design and LFT&E process. Ship designers can benefit by 
gaining insight into the ship's structural response early in the ship design process when it 
is most useful. Also, these simulations can be used to extrapolate shock trial data. This 
could lead to knowledge of the ship's response to realistic threat conditions or the effect 
of configuration changes on the ship's response. This last item is significant since our 
ships have historically undergone extensive modifications, structural and otherwise, 
during their lifetimes. Finally, increased concern over the environment has lead the Navy 
to invest large amounts time and money into preparing for, and evaluating the effect of, 
ship shock trials. Modeling and simulation could help reduce the number and size of 


UNDEX events required during live fire ship shock trials [Ref. 5]. 
E 


B. SCOPE OF RESEARCH 

This paper reviews previous research done at NPS on modeling and simulation of 
the response of a surface ship subjected to an underwater explosion. These previous 
research efforts provided the background for the modeling and simulation of the ship 
shock trials of the USS John Paul Jones (DDG-53) [Ref. 7]. One- and Two-Dimensional 
models are presented as a method of validation for the Three-Dimensional work. It 
further investigates the optimum depth of the three-dimensional fluid volume that must 
be modeled to accurately capture the effect of cavitation on the surface ship structural 
response. Any reduction in the size of the fluid model cuts down on computational 
expense. The 3D model to be considered in this work is a simple rectangular barge-like 
box structure. The effect of varying the fluid mesh depth on selected nodal velocities will 
be compared. The structural and fluid models are constructed using the TrueGrid finite 
element mesh generation program [Ref. 8]. The analysis of the model response is 
conducted using the LS-DYNA/USA (Underwater Shock Analysis) coupled computer 
code [Ref.’s 9 and 10]. Finally, a second generation finite element fluid model for the 
ship shock simulation of the USS JOHN PAUL JONES (DDG 53) will be presented in 
the Appendix. 
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Il. UNDERWATER EXPLOSIONS 


A. UNDERWATER SHOCK PHENOMENA 

When a high explosive, such as TNT or HBX-1 is detonated, the original solid 
material of the explosive is converted into a very high temperature and pressure gas 
within nanoseconds (on the order of 3000°C and 50000 atm.) [Ref. 11]. The pressure 
wave that is formed originates in one section of the explosive and propagates throughout 
the remainder of the explosive. As this pressure wave propagates, it initiates the 


chemical reaction that creates more pressure waves. 






High Pressure Gas Bubble 


Charge Diameter 


Shock Wave Front 


The Figure 1:Gas Bubble and Shock Wave [from Ref. 11] 


The pressure wave velocity steadily increases within the solid explosive until it exceeds 
the speed of sound in the explosive, creating a shock wave. The shock wave propagates 
through the solid at a constant speed and then, with the high temperature and pressure 
behind the shock front, into the surrounding medium [Ref. 11]. At detonation, the 
pressure rise produces a steep fronted discontinuous wave, which decays exponentially 
with time as shown in Fig. 2. The duration of the pressure disturbance lasts only a few 
milliseconds. The shock wave near the charge 1s assumed to propagate at several times 
the speed of sound in water, it then falls rapidly to the acoustic velocity, approximately 


5,000 ft/sec, as it travels radially outward. 


34,000 psi 


3,400 psi 
и 2.200 psi qu 
340 psi 
5 45 50’ 495" 500 


Figure 2. Shock Wave Profiles From a 300 lb. TNT Charge [from Ref. 11] 


Additionally, the pressure profile of the shock wave is proportional to the inverse of the 
distance from the charge, l/d, and the wave profile gradually broadens as it spreads out 
[Ref. 8]. Empirical equations have been determined to define the profile of the shock 
wave. These relations enable calculation of the pressure profile of the shock wave (P(t)), 
the maximum pressure of the wave (Pmax), the shock wave decay constant (0), the bubble 


period (T), and the maximum bubble radius (A max). 
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Other variables in the equations are: 


W = Charge weight (Ibf) 

R = Standoff distance (ft) 

D = Charge depth (ft) 

tı = arrival time of shock wave (msec) 

t= time of interest (msec) 

Ky, Ko, Ks, Ke, A1, A2 = Shock wave parameters 


Through calculation, tt can be determined that Pmax decreases by approximately one-third 


after one decay constant. 


The high-pressure gas that results from the explosion expands outward in a radial 
manner as shown in Fig. 1, and imparts an outward velocity on the surrounding water. 
Initially, the pressure is much greater than the atmospheric and hydrostatic pressure that 
opposes it and is compressive in nature. Subsequent pressure waves or bubble pulses are 
generated by the oscillation of the gas bubble created by the underwater explosion. The 
peak pressure of the first bubble pulse is approximately 10-20% of the shock wave, but is 
of greater duration making the area under both pressure curves similar [Ref. 11]. The 
inertia of the bubble causes it to expand beyond dynamic equilibrium. This means the 
internal pressure is less then the hydrostatic plus atmospheric. The bubble then contracts 
beyond dynamic equilibrium, followed by another expansion. This oscillation sequence 
continues until the energy of the reaction is dissipated or the bubble reaches the free 
surface or impacts the target. Bubble pulse excitation significantly effects the wiping 


response of a ships hull. 


Depending on the charge location, depth of the water and bottom composition, 
other effects are characteristic of an underwater shock. Bottom bounce is the reflection 
of the shock wave off of the bottom of the body of water. If the bottom was perfectly 
rigid a compressive wave originating from an image charge would result. Since no 
bottom composition is near rigid the resulting pressure pulse is much smaller then the 
above assumption predicts. In reasonably deep water, this effect is not an issue for 


surface vessels. 


Free surface reflection is a very important effect. The reflection of a pressure 
wave from a free surface results from the requirement that the pressure above the surface 
be unchanged. and the reflected wave must therefore be one of negative pressure. This 
type of pressure wave is called a rarefaction wave. This rarefaction wave causes bulk 


cavitation. 


B. FLUID-STRUCTURE INTERACTION 

When an object such as a ship is in the vicinity of an underwater explosion, the 
shock pressure pulses produced by the explosion impinge upon the surface of the 
structure. A fluid-structure interaction takes place as the pressure pulse acts upon the 
flexible surface of the structure. This discretized equation of structural motion can be 


expressed by: 
IM IS + [C, 19 + IK, 199 = (Ej (2.6) 


where [Ms] is the mass matrix, [Cs] is the dampening matrix, [K;] ıs the stiffness matrix, 
{x} is the acceleration vector, {x} is the velocity vector, and {x} is the displacement 
vector of the structure and {f} is the external force vector. In the case of a submerged 


structure excited by an acoustic wave, íf) is given by: 


{В = [ОЛА ри} + {р5}) + ths (2.7) 


where [G] is the transformation matrix that relates the surface nodal forces of the fluid 
and structure, [Af] is the diagonal area matrix associated with the fluid elements, 
fp¡)=Incident wave nodal pressure vector, and {ps}=scattered wave nodal pressure vector 


[Ref. 11]. 


IM, ]{Ps 3 + pe[lA, Jtps 3 = pelM, ltüs} (2.8) 


The approximate relation in eqn. 2.8 ıs called the Doubly Asymptotic Approximation 
(DAA) method because it approaches exactness at both low and high frequencies. This 
equation describes a matrix of differential equations in time for the approximation of 


acoustic fluid-structure interaction. The DAA represents the surrounding fluid of the 


structure through the interaction of state variables pertaining only to the structure’s wet 


surface. The terms of the equation are: [Mr] 1s the symmetric fluid mass matrix for the 
wet-surface fluid mesh, (ps) and (ps) are the nodal pressure vector and its first time 
derivative of the scattered wave, c 1s the acoustic velocity of water, [Af] 1s the diagonal 
area matrix associated with the fluid elements, and (ús) is the scattered wave velocity 


vector. [Ref. 10] 


The kinematic compatibility relation can then be applied to relate (us) to the 


structural response, 


[G]" {x} = {u,} + {us} (2.9) 


The “T” superscript indicates the transpose of the matrix. This equation is an expression 
of the constraint that the normal fluid particle velocity must match the normal structural 


velocity on the structure wetted surface. 


Substituting Equation (2.7) into (2.6) and Equation (2.9) into (2.8) results in 
Equations (2.10) and (2.11), 


ІМ, (Xj € [C] o5 + [К Кх = СПА, (р, + {Ps}) (2.10) 


апа 


[M; ] {Ps} + pc[A, J{ps} = рМ, (61 (33 - (41)) (2.11) 


The Underwater Shock Analysis (USA) code solves Eqs. (2.10) and (2.11) 
simultaneously by using a staggered solution procedure that is unconditionally stable 
with respect to the time step used [Ref. 10]. Once this system of equations is solved, 


desired response results such as displacement, velocity, and acceleration can be studied. 


C. CAVITATION 
Two types of cavitation can occur during an UNDEX event. "Local cavitation" 
occurs at the fluid-structure interface and "bulk cavitation" occurs near the free surface 


and can cover a relatively large area. 


1. Local Cavitation 


Taylor flat plate theory, the simplest case of fluid-structure interaction is used to 
illustrate how local cavitation occurs. In this case, an infinite, air-backed plate is acted 


upon by an incident plane shock wave as shown below in Fig. 3. 
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Figure 3. Taylor Plate Subjected to a Plane Wave [Ref. 9] 


Once the shock wave strikes the plate, a reflected shock wave leaves the plate. 


According to Newton s second law of motion as shown in Eq. (2.12), 


ii (2412) 


dt 


where m is the mass of the plate per unit area, u is the velocity of the plate after being 
subjected to the shock wave, P;(t) is the incident wave pressure and P(t) is the reflected, 
or scattered, wave pressure. Define the fluid particle velocities behind the incident and 
reflected shock waves as u;(t) and up(t). The velocity of the plate is then defined by Eq. 
(ІЗ 


u(t) = u,(t)—u,(t) (2.13) 


For the one dimensional plane wave, the wave equation is P = pCu. It follows that the 


incident and reflected shock wave pressures, Eqs. (2.14) and (2.15), are 


P, 2 pCu, (2.14) 
P, = pCu, (2.15) 


where p is the fluid density and C is the acoustic velocity in water. Substituting the 
above pressure Eqs (2.14) and (2.15) into the velocity Eq. (2.15) results in the incident 


shock pressure. The incident shock wave pressure is defined as Eq. (2.16) 


P(t)=P, e? (2.16) 


Solving for the reflected shock pressure yields Eq. (2.17), 


Р,(@)=Р,—рСи=Р „е ° —pCu (2.17) 


where t is the time after the shock wave arrives at the target. Now the equation of 


motion, Eq. (2.12) can be rewritten as Eq. (2.18), 
du E 
ne) +pCu=2P е (2.18) 


which is a first order, linear differential equation. The solution, u(t), of the differential 


equation IS expressed In Ea. (2.107 as 
Bt 
E 


о 0 
u = — _ 
m(1—B) 


with B = pC@/m and t>0. The total pressure that impinges on the plate is defined as Eq. 
(2.20), 


[e ° — ° °] (2.19) 
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As the value of B becomes larger, as in the case of a lightweight plate, the total pressure 


will become negative at a very early time. However, since water cannot support tension, 


negative pressure cannot exist. Therefore, as the water pressure reduces to vapor 


ІМ 


pressure at the surface of the plate, cavitation occurs. At this point, the pressure in front 


of the plate has been cut off and the plate has reached its maximum velocity. [Ref.12] 


A ship's hull can be easily generalized as a Taylor flat plate. Local cavitation is 
likely to occur along the hull where the pressure pulse from the UNDEX impinges with 
sufficient force and the hull plating f value is large enough to make the net pressure 


negative. 
2. Bulk Cavitation and Vertical Kick-off Velocity 


As detailed above, a rarefaction wave, which is tensile in nature, is created when 
the shock wave is reflected from the free surface. Since water cannot sustain a significant 
amount of tension, cavitation will occur when the pressure drops to zero or below. Upon 
cavitation, the water pressure rises to the vapor pressure of water, approximately 0.3 psi. 
This cavitated region created by the rarefaction wave is known as the bulk cavitation 
zone. It consists of an upper and lower boundary and its extent is dependent on the charge 


size, type, and depth [Ref.’s 11 and 12]. 


Figure 4 shows a typical bulk cavitation zone. The cavitation zone is symmetric 
about the y-axis in the figure. The water particles behind the shock wave front at the time 
of cavitation have velocities depending on their location relative to the charge and the 
free surface. Water particles near the free surface, for example, will have a primarily 
vertical velocity at cavitation. As the reflected wave passes, the particles will be acted 


upon by gravity and atmospheric pressure. 


An estimate of the vertical kick-off velocity of a surface ship can be arrived at by 
looking at the fluid particle velocities near the free surface. Since the mass of the ship is 
the same as the water it displaces, the water particle velocities through this volume can be 
calculated to give the ships velocity. The centerline can be calculated and used as the 


average vertical kickoff velocity of the ship section. [Ref. 12] 


The upper cavitation boundary is the set of points where the rarefaction wave 
passes and reduces the absolute pressure to zero or a negative value. The region will 
remain cavitated as long as the pressure remains below the vapor pressure. The total or 
absolute pressure, which determines the upper boundary, is a combination of atmospheric 


2 


pressure, hydrostatic pressure, incident shock wave pressure, and rarefaction wave 


pressure. 


Upper Cavitation Boundary 









Free Surface 


Cavitation Region 


Lower Cavitation Boundary 


y 
Figure 4. Bulk Cavitation Zone [Ref. 12] 


The lower cavitation boundary is determined by equating the decay rate of the 
breaking pressure to the decay rate of the total absolute pressure. The breaking pressure 
is the rarefaction wave pressure that reduces a particular location of a fluid to the point of 


cavitation pressure, Or Zero psi. 


The upper and lower cavitation boundaries can be calculated from Equations 
(2.21) and (2.22), respectively [Ref. 12]. Any point which satisfies F(x,y) and G(x,y) = 0 


determines the bulk cavitation boundary. 
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The variables in Equations (2.21) and (2.22) are: 
x, y = horizontal range and vertical depth of the point 
rı = standoff distance from the charge to the point 
r2 = standoff distance from the image charge to the point 
C = acoustic velocity in the water 
D = charge depth 
8 = decay constant 
y = weight density of water 
Pa = atmospheric pressure 
W = charge weight 
P; = P(t), Equation (2.1) 
9 = Equation (2.3) 
Ki, Aj 7 shock wave parameters 


Figure 5 shows the charge geometry for the above two equations. 


"m Image Charge 





Figure 5. Charge Geometry for Bulk Cavitation Equations [Ref. 12] 
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Appendix A provides a MATLAB m-file that calculates and plots the bulk 
cavitation zone for a user supplied charge weight and depth by solving Equations (2.21) 


and (2.22). 
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III. PREVIOUS NPS UNDEX RESEARCH 


Significant work in the area of surface ship modeling and simulation has been 
accomplished at the Naval Postgraduate School. The work of Santiago [Ref. 14] and 
Wood [Ref. 13] provided the background experience and modeling techniques that led to 
promising results in the modeling and simulation of the ship shock trial of USS John Paul 
Jones (DDG-53) as presented in Ref. 7. 

Subsequent developments in the modeling and simulation codes utilized in Ref.”s 
13 and 14 have greatly simplified these tasks. To illustrate the streamlining of the 
modeling and simulation process consider the number of different computer codes 
utilized in the previous work. In Ref. 14, Santiago used 2 modeling codes, 6 analysis 
codes and 4 post-processors. In Ref. 13, Wood used 2 modeling codes, 3 analysis codes 
and 4 post-processors. This study uses 1 modeling code, 2 analysis codes and 2 post- 


processors. 


A. ONE-DIMENSIONAL MODEL 

Santiago in Ref. 14 offered an example of a one-dimensional problem involving 
the interaction between an air-backed plate and a cavitating fluid. This analysis utilized 
the USA-CFA (Underwater Shock Analysis-Cavitating Fluid Analyzer) [Ref. 15] code to 
solve this problem. The excitation is from a plane exponential wave impinging upon the 
plate. H.H. Bleich and I.S. Sandler, describe an exact solution to this problem [Ref. 16]. 
Ref. 14 produced identical results as those shown by Bleich and Sandler. 

The method utilized by Santiago has been greatly simplified due to subsequent 
development of the USA and LS-DYNA codes. The one-dimensional model was 
constructed using the TrueGrid mesh generation code. All nodes in the model are 
constrained allowing only freedom in the X (vertical) direction. The model represents an 
air-backed plate at the free surface of a cavitating fluid. The structural plate is modeled 
by a single 1.5-in square shell element. The structure was constructed of 1-in steel plate 
modeled by LS-Dyna material Type 1:Mat Elastic (weight density of 0.284-Ibf/in?, 
Young's Modulus of 30E6-psi and Poison's Ratio of 0.3) [Ref. 9]. The fluid column is 


I 


modeled by 100 hexahedral solid elements and contains 404 nodes. The fluid is modeled 
by LS-Dyna material Type 90:Mat_Acoustic (weight density of 0.0361-Ibf/in’, acoustic 
speed of 4760-ft/sec) [Ref. 9]. The fluid volume was constructed using the TrueGrid’s 
BLUDE command. This command extrudes a set of polygons through the mesh lines of 
a user-defined block (grid) [Ref. 8]. The USA and LS-DYNA input decks used for this 
model are contained in Appendix B. Figure 6 shows the model geometry. The DAAI 


boundary for this model 1s defined at the bottom of the fluid column. 
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Figure 6: 1D Plate and Cavitating Fluid Model 
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B. TWO-DIMENSIONAL MODEL 

Wood in Ref. 13 investigated the effect of cavitation on a ship-like box structure 
subjected to an UNDEX event. The study looked at two- and three-dimensional models 
and compares the effect of modeling with and without a coupled fluid volume. 


Additionally, Ref. 13 investigates the effects of togeling the LS-Dyna material Type 


18 


90:Mat Acoustic cavitation flag on and off and Rayleigh damping. The LS-Dyna and 
USA input decks for the present study are based largely on those developed in Ref. 13. 

Wood in Ref. 13 presents a two-dimensional model to perform the initial analysis 
and verify the shock wave behavior in the fluid mesh. The two-dimensional model 
consists of a midships section of a ship-like box structure and a surrounding fluid 
volume. A similar model is presented here to illustrate the shock wave propagation in the 
models. 

The two-dimensional model was constructed using the TrueGrid mesh generation 
code. The structural model represents a cross-sectional area of a barge. It is 24-in wide, 
24-in high and 6-in deep. It contains two symmetric boundary conditions on the XY- 
planes. The shell plating was constructed of l4-in steel plate modeled by LS-Dyna 
material Type 1:Mat Elastic (weight density of 0.284-Ibf/in?, Young's Modulus of 30E6- 
psi and Poison's Ratio of 0.3) [Ref. 9]. The finite element mesh of the barge structure 
contains 38 nodes and 18 shell elements. The fluid volume is 80-in wide, 120-in high 
and 6-in deep. The fluid is modeled by LS-Dyna material Type 90:Mat_ Acoustic (weight 
density of 0.037-Ibf/in’, acoustic speed of 4929.6-ft/sec) [Ref. 9]. The fluid volume was 
constructed using the TrueGrid's BLUDE command. The finite element mesh of the 
fluid volume contains 616 nodes and 270 hexahedral solid elements. Figure 1 shows the 
fluid volume finite element model. 

Of particular importance in 2D models utilizing the USA code is the numbering 
of nodes along the DAAI boundary. This boundary is located around the circumference 
of the fluid mesh. For example, the DAAI boundaries are the faces parallel to the Z axis 
in figure 10, excluding the free surface. Referring to figure 8 and imagining the 2D 
cross-section as a cylinder section, the first two nodes of these DAAI face elements 
should progress around the circumference not along the Z axis. The NSORDR command 
in the USA FLUMAS input deck allows the user to cyclically reorder these nodes to 
achieve this requirement [Ref. 17]. See the 2D USA FLUMAS input deck in Appendix 
B. 
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Figure 7: 2-D Fluid Model with Barge Cross-section 
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Figure 8: 2D Faceset Node Numbering 
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IV. MODELING AND SIMULATION 


Modeling and simulation involves model construction and pre-processing, 
analysis and solution, and post-processing programs. A flow chart of the model building 


and testing procedure is shown below in Figure 9. 


Model Construction | 
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Figure 9: Flow Chart. Model Construction and Simulation 


A. MODEL CONSTRUCTION AND PRE-PROCESSING 
1. 3D Structural Model 


The three-dimensional model was constructed using the TrueGrid mesh 
generation code[Ref. 8]. The structural barge model is 168-in long. 24-in wide. and 24-in 
deep. It contains two structural athwartships bulkheads. The model was loaded with 
three lumped masses (0.664-lbf s /in^) evenly spaced on the centerline to give the model 
a 12-in waterline. The shell plating was constructed of %-in steel plate modeled by LS- 
DYNA material Type 1:Mat_Elastic (weight density of 0.284-Ibf/in’, Young’s Modulus 
of 30E6-psi and Poison’s Ratio of 0.3)[Ref. 9]. The finite element mesh of the barge 
structure contains 443 nodes and 452 shell elements. Figure 10 shows the barge finite 


element model. 


Athwartships 
Bulkheads 





Figure 10: Barge Finite Element Model 


2. 3D Fluid Modeling 


The fluid is modeled by LS-DYNA material Type 90:Mat_Acoustic (weight 
density of 0.037-Ibf/in?, acoustic speed of 4929.6-ft/sec). The fluid volume was 
constructed using the TrueGrid's BLUDE command. This command extrudes a set of 
polygons through the mesh lines of a user-defined block (grid). Figure 11 shows an 
example of one of the fluid volume finite element models used in this study. The finite 
element mesh of this fluid volume contains 18615 nodes and 16272 hexahedral solid 


elements. 


The nodal spacing adjacent to the structural model is critical for the stability of the 
USA analysis. The nodal distance normal to the structural mesh limits the size of the first 
layer of fluid elements [Ref. 15]. This parameter is determined by the following 


equation: 


20025 


Pst, 
Where p = density of water, D = thickness of the fluid element in the direction normal to 
the wetted surface of the structure, p, — density of the submerged structure, and t, = 
thickness of the submerged structure. For the barge model the critical element thickness, 
D, is 4.8 inches. The first layer of fluid mesh for these models was set to 4-in. The bulk 
fluid elements were generated using the RES and AS commands. These commands space 
the nodes of a mesh so that the ratios of the distances between them will be constant [Ref. 
8]. This means the fluid elements gradually increase in size as you move away from the 
structure. These commands were used to generate a consistent mesh quality at a given 
distance from the structure. The length and width of each fluid volume are identical. The 
depth of each model is varied Ітоп 324-іп in the Base model to 42-in in the Third model. 
Table 1 summarizes the specifics of each model. Figure 12 shows the truncation of the 


fluid volume for the Base, Extended, Full and Half models. 
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Table 1:Fluid Volume Models Specifics 
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This thesis also presents a second-generation coupled structural-fluid model of the 


USS JOHN PAUL JONES (DDGS53). A detailed description of the DDG53 model is 


contained in Appendix D. 
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Figure 11: Example of Fluid Volume Finite Element Model 
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Figure 12: Fluid Volume Truncation 


B. ANALYSIS AND SOLUTION 


1. Analysis Program Description 


The finite element model must be translated into LS-DYNA keyword format 
using the TrueGrid Output translator. The USA code performs the majority of the work 
(formulation of the fluid-structure interaction matrices) and LS-DYNA is used to perform 
the time integration solution for the structure. LS-DYNA is a non-linear three- 
dimensional structural analysis code [Ref. 9]. The USA code consists of three main 


modules: FLUMAS, AUGMAT, and TIMINT [Ref. 10]. 


FLUMAS is the first USA module required to be run. FLUMAS generates the 


fluid mass matrix for the submerged portion of the structure. The fluid mesh data, as well 
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as the transformation coefficients that relate both the structural and fluid degrees of 
freedom on the wetted surface are generated, including the nodal weights for the fluid 
element pressure forces and the direction cosines for the normal pressure force. The fluid 


area matrix 1s diagonal and the fluid mass matrix is fully symmetric.[Ref. 10] 


AUGMAT is the second module to be run. This module takes the data generated 
by the FLUMAS and initial LS-DYNA runs to construct specific constants and arrays 
utilized in the staggered solution procedure for the actual transient response analysis 
[Ref. 14]. The AUGMAT module organizes this information so that the integrator need 
refer to only one permanent data file for input. The augmented interaction equations are 


formed from Equations (2.10) and (2.11). [Ref. 10] 


TIMINT performs the direct numerical time integration and also handles the 
computation of the UNDEX parameters, such as the shock wave pressure profile. The 
structural and fluid response equations are solved separately at each time step through the 
extrapolation of the coupling terms for the two systems. LS-DYNA is used to solve the 
structural equations and the TIMINT run solves the fluid equations. A result of using the 
staggered solution procedure mentioned previously is that LS-DYNA and TIMINT can 
each have a different time step assigned. However, the general practice is to use the same 
time step in both computations. Despite using an unconditionally stable solution scheme, 
the TIMINT time step must be set small enough to accurately capture the fluid system 
response. Additionally, it should be mentioned that LS-DYNA uses a central-difference 
integration method, which is conditionally stable. The LS-DYNA time step must be less 
than or equal to the critical time step for the structural finite element mesh or numerical 


instability will result. [Ref.’s 9 and 10] 


Appendix B provides example input decks for each of the three USA modules as 
well as an example LS-DYNA Keyword input deck. 


2. 1D Model 


Figure 13 shows the pressure time history of the fluid element immediately below 


the plate for the 1D case. The plots show the Cavitation flag on and flag off cases. These 
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results show excellent correspondence with those in Ref. 14 and Ref. 16. This confirms 


the validity of the abbreviated modeling and simulation techniques. 
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Figure 13: Pressure and Velocity Response of 1D Model 


3. 2D and 3D Test Description 


The charge placement, designated SHOT-1, used in the shock simulation run for 


the 2D and 3D models in this study was a port side beam on shot. A 60 Ib. TNT charge 
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was selected for the simulations. SHOT-1 placed the charge offset from the side of the 
barge at a horizontal standoff distance of 10-ft (120-in) and a depth of 25-ft (300-in) from 
the free surface. The coordinate location of the charge was (84,-132,-288-in). Figure 14 
shows this attack geometry and Table 2 shows a summary of the UNDEX parameters of 
the explosion. As listed in Table 2, the bubble period is over 500-msec so there is no 
bubble interaction in the 50 msec run-time of this study. Figure 15 shows the theoretical 
pressure profile at the barge for this shot geometry. The times listed in Figure 15 are 
taken from the instant of detonation. The bulk cavitation zone was computed and is 
shown in Figure 16. The theoretical depth of the cavitation under the barge structure is 
84-in. Since the mass of the surface ship 1s the same as the water it displaces, an estimate 
of the vertical kick-off velocity of the ship can be found based on the water particle 
velocities [Ref. 12]. For the geometry of SHOT-1 the theoretical kick-off velocity was 


calculated as 20.7-ft/sec. 
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Figure 14: SHOT-1 Geometry 


Table 2: UNDEX Parameters 


Pm | 2409 psi | 
e | 0322msec 
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Figure 16: Cavitation Zone for SHOT-1 


4. 2D Model Results 

The pressure time histories of three fluid elements are plotted in figure 17. These 
elements represent the fluid pressure in the top, middle and bottom of the coupled fluid 
volume. The magnitude and decay of these plots matches what is expected from the 
geometry of the UNDEX attack. The second peaks seen are due to the reflection of the 
incident shock wave off the structure. The cut off of the pressure in the top and middle 
elements represents cavitation. The bottom element does not experience this cut off 
indicating no cavitation at this depth. This confirms the proper behavior of the model 
since the theoretical cavitation depth was 84-in. Figure 18 shows the incident shock 


wave in the 2D fluid mesh. 
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Figure 17: 2D Pressure Response 
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Figure 18: 2D Incident Shock Wave 
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(с: POST-PROCESSING 

The solution data is output into two formats for analysis: binary and ASCII. The 
binary data files created by the LS-DYNA/USA runs contain the model’s finite element 
response information. LS-POST [Ref 18] was used for three-dimensional response 
visualization. It provides animation and image generation features. This allows the 
animation of the shock wave propagation through the fluid to the structure. Additionally, 
LS-POST has the capability of extracting ASCII solution data, plotting time histories of 
fluid and structural responses (pressure, stress, velocity, acceleration, etc.) and writing it 
to a separate ASCII file for later evaluation. The Time history plots were made by 


exporting the LS-POST data to Excel. 
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V. 3D SIMULATION RESULTS 


This study examined the effect of decreasing the fluid mesh depth on the nodal 
response of the barge structure. The five fluid models studied were explained above and 
are simply labeled Base, Extended, Full, Half and Third. The nodes analyzed are listed in 
Table 3 and shown in Figure 19. These nodes were chosen based on their physical 


location in the structure. 


Table 3: Response Nodes Coordinates 


„0, 0, (84,0,0) 


x,y,z Coordinate | (0,0,12) (50,0,12) 
(inches) 
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Figure 19: Nodes Analyzed 
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The vertical response of the selected nodes are calculated and compared for the 
five fluid volume models. The Base model is used as a reference for the comparisons. 
The time history plots of the Base and Third models for the selected nodes are offered as 
an example of the response data. Although they were developed to compare measured 
verse calculated data, the Russel and updated Geers error measurement criteria are 
adapted in this study to quantify the effect of truncating the fluid volume depth [Ref.”s 
19,20 and 21]. The formulation of these error factors is shown below. 
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Updated Geers Error Formulation [Ref. 21] 
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Russel Error Formulation [Ref. 19] 


The tabulated results of the magnitude, phase and comprehensive error factors are 


presented. The Russel error factors are designated RM, RP and RC. Reference [19] 


suggests if the magnitude, phase and comprehensive Russel error factors are less then 0.2 
the correlation between the data is acceptable. However, since this is a sensitivity 
analysis of simulation models a factor of 0.1 is used. The updated Geers error factors are 


designated GM, GP and GC. 
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Figure 20: Resulting Pressure Profile at Midships of Barge 


A. INCIDENT SHOCK WAVE AND BULK CAVITATION FORMATION 


Figure 19 shows the theoretical and simulation results pressure-time history of 
fluid (brick) element 2239 of the Base model. This element is located on the port side of 
the structure at coordinate (84,-16,-4). This pressure profile is typical of all four models 
and compares well to the theoretical profile. As seen in Figure 19, the pressure rises 
quickly to a peak, dips and then peaks again at approximately 0.35-msec. This second 
peak is due to reflection of the incident pressure wave off the structural model. The 
pressure drops off to zero at 0.6-msec representing cavitation. This pressure-time 
response shows reasonable correlation with the 2D model presented above. Figures 21- 


25 show the incident shock (fig. 21), initial cavitation bubble formation (fig. 22-23) and 


on 


complete envelopment of the barge by bulk cavitation (fig 24-25). Figures 22 and 23 are 


— 
er 


at time 0.4-msec and figures 24 and 25 are at time 1.09-msec. These events correspond 


well with the pressure-time history of figure 19. 
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Figure 23: Top View of Initial Bulk Cavitation Zone 
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Е igure 24: Barge Enveloped by Bulk Cavitation 





Figure 25: Top View of Barge Enveloped by Bulk Cavitation 


B. VERTICAL VELOCITY RESPONSE 


Figure 26 shows the vertical bodilv response of the barge model to the UNDEX 
attack. The results plotted represent the vertical velocity of the barge structure from the 
time of impact of the incident UNDEX shock wave out to 50-msec. The response is as 
expected from the phvsics of the attack. The incident shock wave 1mpacts the structure 
with extremely high pressure and forces the barge rapidlv upward. This initial velocitv of 
approximately 2! ft/sec corresponds well to the theoretical kickoff velocity of 20.7 ft/sec. 
The velocity continues to increase until approximately 1-msec, the onset of bulk 
cavitation. At 4-msec the velocity of the structure hits the peak value. approximately 23- 
ft/sec. Тһе barge velocity then drops off and becomes negative at 45-msec. The 
oscillations in the time history are due to the dynamic structural response. such as 
bending and torsional modes, of the barge. (For an insight into the effect of the coupled 
fluid volume, see Appendix E for a comparison of coupled verse uncoupled model 


responses). 
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Figure 26: Vertical Response of Barge Model 
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Figure 27: Vertical Response of Node 25 for Base and Half Models 
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The vertical response of node 25 for the Base and Third models are plotted in 
Figure 27. This node 15 located at (0,0,12) in the center of the forward bulkhead at the 
extreme forward end of the barge. Due to its location, the beam bending of the barge 
structure as well as the plate buckling of the bulkhead influence the vertical response of 
this node. Initially, the velocity rises sharply due to the impact of the incident shock 
wave. There is then a sudden drop in the velocity due to the bending response of the 
structure. This bending response combines with the overall bodily response to cause a 
high peak velocity of 37 ft/sec. The effect of reducing the fluid volume depth is 
immediately apparent at approximately 12-msec. At this time the Third model separates 
from the Base response curve. As Table 4 shows, the phase, magnitude and 
comprehensive error factors all compare very well for all models analyzed except the 
Third model. The response is out of tolerance for the Third model based on the Geers 
magnitude and comprehensive error factors. 


Table 4: Russel’s and Geer's Error Factor for Node 25 
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The vertical response of node 159 is plotted in Figure 28. This node is located at 
(50,0,12) in the center of the forward athwartships bulkhead. Table 5 presents the error 
factors for this node. Both Russel and Geers error factors are above tolerance at this 


node. 
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Table 5: Russel’s and Geer's Error Factor for Node 159 





[Extensa] Fan | Har ти” 
aS Aa. a. © 
RP 008 

WoWWedGees| | | T 
018 































Time (msec) 


— Base — Third 
40 - r 
| 
| | | | | 
ih | dc 
қ A л A к | та IE | \ | : 
| 1 M | | || т | ] AA 
| | | ' y | Ç YA А Ë ATA", | ja | А, | i 
| | | | | | | № 
E | | ! — m. on | 
| | | | | 
zip ue ы A | T a |с 2 
0 5 10 15 20 25 30 35 40 45 


Figure 28: Vertical Response of Node 159 for Base and Half Models 
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Figure 29: Vertical Response of node 285 for Base and Half Models 
The vertical response of node 285 is plotted in Figure 29 and the results of the 
error analysis in Table 6. This node 15 located at (118,12,12), on the starboard side of the 
barge at the intersection of the aft athwartships and side bulkheads. Due to its location, 
this node is not as affected by the bending response of the barge structure as node 25. 
Additionally, the stiffness of the location prevents excessive buckling. The response 
follows that of the overall bodily response, Figure 8, with a high frequency addition. The 


overall frequency response of the models compares well throughout the run. 


Table 6: Russel's and Geer's Error Factor for Node 285 


Extended] Fun] Far] Third 
Rasse | 
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The vertical response of node 213 is plotted in Figure 30. This node is located at 
(84,0,0), оп the centerline of the barge in the middle of the floor. This node is located at 
very flexible point in the structure and at the location of a lumped mass. These two 
factors damped the response of the structure at this node. The high frequency responses 
observed in the previous nodes are absent here due to increased flexibility at this location. 
The negative sign on the magnitude error factors indicate that the magnitude of these 


models was greater then that of the Base model. 


Table 7: Russel's and Geer's Error Factor for Node 213 
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Figure 30: Vertical response of Node 213 for Base and Half Models 


Figure 31 shows the vertical response of node 398. This node is located at 
(142,12,12) on the starboard side bulkhead. This node is located at a very flexible point 
in the structure and is subjected to plate buckling. Again, the velocity rises sharply and 
then drops off due to the bending response of the structure. As seen in table 8 the error 


factors are well under the designated criteria. 


Table 8: Russel's and Geer's Error Factor for Node 398 
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Figure 31: Vertical Response of node 398 for Base and Half Models 
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C. ATHWARTSHIPS AND LONGITUDINAL RESPONSE 

The major excitation on the barge model due to the geometry of this UNDEX 
attack is vertical. For completeness the longitudinal and athwartships velocity response 
of several nodes is included. It should be noted that the maximum velocity in these 


directions is less then 20% of the vertical. 


The longitudinal response of node 285 is plotted in Figure 32 and the results of 
the error analysis in Table 9. This node is located at (118,12,12), on the starboard side of 
the barge at the intersection of the aft athwartships and side bulkheads. The stiffness of 
the location prevents excessive buckling. As seen in figure 32 there is poor magnitude 
and phase correlation between the Base and Third models. Tables 9 and 10 show there is 


no correlation between all the models for the longitudinal response. 
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Figure 32: Longitudinal Response of node 285 for Base and Half Models 
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Table 9: Node 285 Longitudinal Velocity 
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Table 10: Node 213 Longitudinal Velocity 
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Figure 33 shows the athwartships velocity response of node 285. The Base and Third 
models separate at approximately 15-msec. There is little magnitude and phase correlation 
between these models. As shown in Tables 11 and 12 all the models experience large errors for 


the magnitude, phase and comprehensive error factors. 
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Figure 33: Longitudinal Response of node 285 for Base and Half Models 


Table 11: Node 285 Athwartships Velocity 
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Table 12: Node 213 Athwartships Velocity 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


This thesis reviewed some of the past work done at NPS in surface ship UNDEX 
modeling and simulation and demonstrated the application of updated software to these 
problems. The significant advances in computer codes have greatly simplified the task of 
modeling and simulation. Additionally, this study examined the effect of decreasing the 
fluid mesh depth on the nodal response of a simple barge structure subjected to an 
underwater explosion. The time history responses of selected nodes are presented. The 
Russel and updated Geers error measure criteria are used to quantify the affect of 
truncating the fluid volume depth. For the vertical response, the Third model was above 
the 0.1 tolerance for all of the nodes analyzed except node 213. In all cases, the 
Extended, Full and Half models had error factors of the same magnitude while the Third 
model was significantly larger. The largest difference in the Russel and Geers error 
factors are seen in the phase response. This is attributed to the fact that the Russel phase 
error factor is more sensitive to changes in small and large differences in phasing, and 
less sensitive in the middle range, compared to the Geers phase error factor [Ref. 19]. 
The longitudinal and athwartships velocity responses showed little or no correlation 
between the models. However, since the primary excitation in an UNDEX event is 
vertical this does not significantly affect the validity of this simulation method. The 
simulations conducted in this study demonstrate that the fluid mesh may be modeled to 
approximately one-half of the theoretical cavitation depth and still provide an acceptable 
vertical response simulation. 

Recommendations for areas of additional study: 

je Work on capturing the full gravitational effects on the barge later ın the 

response. It is suggested that a pressure load can be applied to the 
structure to more accurately simulate the effect of gravity. 

— Work is needed to more accurately capture the longitudinal and 

athwartships response of the structure. Further refinement of the 


structural finite element model may help in this area. 
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Much consideration 1s being given in the Navy to advanced hull forms. 
Utilize this technique to compare the response of various hull forms to 
the incident UNDEX shock wave. Such hull forms as SLICE, wave- 
piercing catamaran, and SWATH could be compared to a traditional 


displacement hull. 


APPENDIX A. BULK CAVITATION PROGRAM 


The following program code calculates the bulk cavitation zone by solving 


Equations (2.21) and (2.22). 


Si) Peo pr Malone 


$Bulk cavitation program 


clear 
Pa = 14.7; 
пие 5000; 


format short g 

ашер( "1 - НВХ-1, 2 - TNT, 3 = Репго ген) 

displi !) 

explosive = input(' Tnpuerthe number Corresponding Commn type or 
explosive: '); 

disp(' г) 

weight = input (' Input charge weight (lbs): '); 
Epuspt' ^") 

D= inout (' Input charge depth(ft): '); 

ISE = 

R = input (' Input standoff distance in Eeet: '); 
Emep(' ') 

Gasp (' ©) 


$ Explosive Paramenters 


disp(' ') 

W-weight; 

if explosive == 
К1=22347.6; А1=1.144; 
K2.056,.A2- 2.2077 
ВА. AO, 
К4=3086.5; А4=2.039; 
К5=4. 761; 
K6=14.14; 

elseif explosive == 
5І-22505; А1-1.18; 
K2=.0539; A2=-. 195; 
K3-21.798; A3-.99; 
К4-3034.9; А4-2.155; 
K5=4.268; 
Ка- 12.67; 

else explosive == 
Kl=24589; Al=1.194; 
к- ше А. 
K3=1.674; A3=.903; 
K4=3135.2; -A4=2. 094; 
K5=4.339; 
K6=12.88; 

end 


gamma = 0.037037; 


o 


¿theta = decay constant 


$x = horizontal distance 

$y = vertical distance 

$r1,R = standoff distance from charge to point 

$r2 = standoff distance from image charge to point 

%$Pi = incident shock wave pressure at tc 

upper = []; $Create matrix to store upper boundary data 
lower = ||; Create matrix to store lower boundary data 


$Calculate Upper Boundary 
О-о TOO 
Tor 540: 227200 
rl Soup DL E 
r2 sore (DY) LEKE 
theta = K2*W> (1/3)? (W732 S tb 22710090: 
ЕЕК (И (зури Altexp(=- (r2 rD епеса)))ғРағ(ааптаут ОТЕ 
Re | Wiel AE АНТ; 


if F «e 0 $Test Гот саутгашоп 
upper = [upper xw 
break,end 
end 
end 


$Calculate lower boundary 
for < = 0: (Upper (lengien (upper), 1)) 
fOr Ty - 0022250 
ll Е Пар холу 
х2 = загі ( (0+у) ^2+х^2); 
Get au zeK zn 1/2 )Z (W173, 22710005 
РЕ уат т АТ tex кои СЕ, 
ОС (Pi (Ci lea) t (]r(((ry2o=(2#D P Dp S)  r2))yrl (ae ee ee 


ӘГІДІСБІ O E ADA A E A ll ODE ee 
(A1/r2)* (P1+Pa+ (gamma*y*12)); 


if G >= 0 %Test for cavitation 
lower = [lower; x v]; 
break,end 
end 
end 


$Find tangent point 

pague tandousperi:-, 2»eloweri.2))57 
Upper—uppern «(1 lengren (tan) ) у; 
lower=lower((l:length(tanu)),:); 


$Plot cavitation boundary 
figure(1) 

max" x—_mas (Upper (2,107 
max yemax(lower(:,2)); 
maxx-num2str(max x); 
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maxy=num2str (max у); 
wt = num2str (W); 
depth = num2str(D); 
plot (upper l), -upper i lower ЕШ у, Ея 
axis([10,1200, MOTOT 
type=explosive; 
if type ==1 
title(["Cavitation Zone for a "Twt, IDEI 1 етае аа ре оо е 
әсерін, Feet): 
elseif type==2 

сісіе Cavitation one ESE Aa ',wt,' lb INT Charge at a Depth of 
г Черен,” test 10), 
else 

title(['Cavitation Zone for a ',wt,' lb Pentolite Charge at a Depth of 
Weacoth,’ fee I): 


end 
label (['feet; where max cavitation width 15- "/пахх,! ft']) 
ylabel(['feet; where max cavitation depth is= ',maxy,' ft']) 
w=W; 


¿Shock Wave equations 
w=W; 

mM sqrt(i(D-y)^24x^2); 
EG SOLE ( (Dty) “2+x 72); 


пера = К2*\^ (1/3) “(^^ 3). В)^А2: 
Enz = КІ (М (ООУ ДЕ; 


ВЕ. (w^. 339/R) ^AS; V *$psi-sec/in^2 
E K1*w^.333*(w^.333/R)^A4; $1b-in/in^3 
ЕЕ 50^. 3353/00433) ^05/6); вес 


ENSE = K6*w^.333/(D*t33)^.333; bft 
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APPENDIX B. USA/LS-DYNA INPUT DECKS 


This section of this appendix provides example LS-DYNA and USA input decks for 
each of the three types of models presented: 1D, 2D and 3D. Reference 18 provides 
information concerning the various input deck variables. The 1D LS-DYNA input is shown 
first and is completely annotated with the input card descriptions. Use this as a reference 


when reviewing the other input decks. 


1D LS-DYNA input deck 


S This input deck contains a complete description of the input cards 
S which are used in the LS-DYNA/USA UNDEX simulations of this study. 
S Reference this deck as you review the input deck for the other models. 
S 

* KEYWORD 

ЕРІТЕ 

One D model 

Setops jebwitime) 

*CONTROL TERMINATION 

15,0, 0,0, 0 

SSets structural time step size 

ONTROL TIMESTEP 

SInitial time step size,0,0,0,load curve ID that limits max time step,0 
Mor —6,0.9,0,0.05.0.0,1,0 

m nro] parrallel processing 

*CONTROL PARALLEL 

9,1 

SDefine load curve to limit max timestep 

*DEFINE CURVE 

Y 

SEES d 0E-5 

ПИ 5.1.0E-5 

SDefines nodes to be stored in NODEOUT file 

*DATABASE HISTORY NODE 

m7, 23,4 

SDefines how often NODEOUT file is written to. 

*DATABASE NODOUT 

2200001 

SDefines how often D3PLOT file is written to. 

*DATABASE BINARY D3PLOT 

0.00001 

SDefines how often D3THDT file is written to. 

DATABASE BINARY DSTHDT 

0003 

SDefines what data is written to Binary databases. 

e TABASE EXTENT BINARY 

ENS 35. 1,1,1;1;,1 

о, 0,0, 0 

SDefines a surface for coupling with USA code. 

*BOUNDARY USA SURFACE 
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SSegment ID, Wet surface flag, SOR 
O 

SDefines explosion properties 
*INITIAL DETONATION 


Pe@oustic boundary, Coordinate location of detonation, detonation time... 


SPeak pressure at ref node, decay constant, coordinate location 
Sof ref node, ref node ID 

IN ADA O, OOO. O 

MU. 0009968 I. 504050, 0.0.5 

S 

S$ NODES 

5 

NODE 

1920000000098 700,07,.000000000E£00,09700 00000 0 0E E007 5, 1 

20 209. 00 000 00E200 O UOOOOOOOOE OO ii 


. €* 9 o o е o o е 


202% 750000000, 1. 5000000077125 00 0 0 0. 0 0 

I04. 50:000000,0.00000000087007,17 509020092070 

*PART 

po acoustic) 

Бекет Tp Section ID, Ма lL TD 

опа зо 

SECTION SOLID 

SDefines properties of brick elements 

90,8 

SDefines properties of acoustic elements 

*MAT ACOUSTIC 

SMat 1 ID, mass density, sound speed, damping factor, cavitaion Е аа a 
Satmospheric pressure, gravitational constant 

ии o 155=05 5. 712E+04 0. 500 1200 12.7, 556.1 

SCoordinates of free surface, directional cosine of free surface 
me TOOT O OO OO 

SDefines element number, part number, and nodes in element 
*ELEMENT SOLID 

mun SDN], 3,4,2 

СООО ШОЕ а. 


оо е ө о o o € ө 


НЕО 4015102 4035374032 3572308 и 

5 

эватс өпе - plate at top of fluid column 
* PART 

Plate 

SPart ID, Section ID, Mat'l ID 
И pue 

SDefines properties of shell elements 
*SECTION SHELL 

SSsection ID, Element form, Shear Тастап 
ооо 

$Shell thickness at node 

о оо 

SDefines element number, part number, and nodes in element. 
"ELEMENT SHELL 

РТО 
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SDefines mat'l properties 

*MAT ELASTIC 

SMat'l ID, Mass density, Youngs modulus, 
1.8 32 9Е-04 „За 000E5:0 О-о ЕЕ 

SDefine the nodes to be constrained. 
SAllows x tanslation only. 

~ BOUNDARY ОРУ ЕТ 
eee) eee ee ee 
*SET NODE LIST 

1 

ос 
СЕО ОЗ 


Poisson's ratio 


401,402,403,404 

5 

SDefines a faceset used for DAA boundary 

[ ОЕТ SEGMENT 

1 

9405. 402 501,0: 0005700. 0. 0080050 DOO DIE-O 
"END 


USA input decks for 1D model. 


FLUMAS DATA FOR 1D Model 


flunam geonam strnam daanam S FLUNAM GEONAM GRDNAM DAANAM 
ШЕГЕ Т = PRIGMT PRITRN PRTAME CALCAM 
ШЕ Е E $ EIGMAF TWODIM HAFMOD QUAMOD 
БИЕ ПТ Е S PCHCDS NASTAM STOMAS STOINV 
ES F T S FRWTFL FRWTGE FRWTGR FRESUR 
БП Е Е S RENUMB STOGMT ROTGEO ROTQUA 
PEE ЕР S PRICOE STRMAS SPHERE OT- iM 
ПИЕТЕ Е 5 OCTMOD CAVELU C ERWTEVETNTCAV 
ше Е 5 ВОТКЕЕ МАЗКЕЕ ВМИНР 

0 404 0 1 5 МӘТЕС "NSTRE NEN NGENF 

Ш 0 0 S NBRA Nee NCAV 

E u5E-D5 5.712E-T04 S RIO CHE 

10 SENVEE 

ши 0. 0. ә РЕРІН СХЕЗ ӘКЕСІН 

13.7 386.088 5 РАТМ GRAVAC 

0 S NSRADE 

0 $ NSORDR 

AUGMAT DATA FOR 1D MODEL 

strnam flunam geonam prenam $ STRNAM FLUNAM GEONAM PRENAM 
DEL PF 5 FRWTGE FRWTST FRWTFL LUMPFM 
LET FT 5 FLUSKY DAAFRM SYMCON DOFTAB 

Е Е ЕЕ з РЕТОМТ РКТТЕЦ РЕВЗШЕНВЕВАЦО 

EFF F $ MODTRN STRLCL INTWAT CFAPRE 

110 $ МТУРОА 

404 1212 3 3 S NSTR NSER NFRE NETR 

I SEK ERES 

PAENE T T + NDICOSAISTART JSTOET JIN® 
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TIMINT DATA FOR 1D MODEL 


prenam posnam S PRENAM POSNAM 
resnam S$ RESNAM WRTNAM 
EEE E S REFSEC FLUMEM 
ЕВЕ Е S TNCSTE CENINT 
840 S NTINT NCHGAL 
Пири -058=6 SS STRIIMZBELTIM 
ТЕ Е Е S EXPWAV SPLINE 
ERDF F S HYPERB EXPLOS 
F P P E S BUBPUL SHKBUB 
1 S NCHARG 
Ds $ HYDPRE 
00.0.0005 Е pis 
АЮ ОООО БЭХ SY 
201 S JPHIST 
To. $ PNORM DETIM 
OES SD ES 
2 $ СНСТУР 
ШЕРІ 795.0 795.0 $ МЕТСНТ 
2000 2000 S NSAVER NRESET 
ЖОО O0 S LOCBEG LOCRES 
EF FE 5 FORWRT STBDA2 
2] 050 070-0 5 ХУ YV 
E SNBISPLA 

2D model LS-DYNA input deck 
ЭКЕТЛОКТ) 
КЕТПЕ 


2-D BOX MODEL WITH FLUID 1/Z-AXIS 
*CONTROL_ TERMINATION 

OPE 0120. 050 

*CONTROL TIMESTEP 

QROOUOT 0909, 02020970, 1710 
*CONTROL OUTPUT 

150, 17 0T. 

*CONTROL PARALLEL 

a d 

*DEFINE CURVE 

1 

UNES 

0.10,0.00001 

*DATABASE HISTORY NODE 

Та PEE 

*DATABASE NODOUT 

0.00001 

*DATABASE BINARY D3PLOT 

0.0001 

*DATABASE BINARY D3THDT 

0.001 

*DATABASE EXTENT 
т 


ВТМАКУ 


— 
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PWACAV ITERAT 
BUOYAN 


VARLIN PACKET 
DOUBDE VELINP 


ZC 
52 


SLANT CHGDEP 


LOCWRT NSTART 
ASCWRT 
ZV 


ООО 

*BOUNDARY USA SURFACE 

Да 

*INITIAL DETONATION 
-1,-288.0,-132.0,6.0,0.0 
2466.0,0.000321,-4.0,-16.0,6.0,448 
$*LOAD SEGMENT SET 

S D maet 

*NODE 

1,0.000000000E-00, -12.0000000,0.000000000E*00, 0,0 
2,0.000000000Е%00,... 


о о ° о о 


e277 -108 .000000, = 401300000007 0500300000006 
Са 1.08..000000,—40.,.0000000 6. 00000008) 0.0 


*PART 
solid 
90 2 90 
BoBC TION SSOL LD 
2 8 


*ELEMENT SOLID 
1590, 39,40. 4p 42; 1, 3,052 


DENDO0,589,585,588, 590, 020819253 6267028 


* PART 
shell 
1 1 1 
*SECTION SHELL 
i 2 0.8333 
0375 0.25 025 0.25 


“ELEMENT SHELL 
ШІП І,3,4;2 

5 

ЖҮП, 3,5,6,4 


oe о o o 


I: 35, 36, 39, 37 
*MAT ELASTIC 
jo a 3e+07 0.3 0 0 


HAT” ACOUSTIC 

ВВ с. ЭЗВ9Е- 4, 59155.24, 10, 14.7,-96.1 

ша Ol. OO. 

BOOUNDARY SPC SET 

Sethis Command is a G r O RI aput a contraint boundary condition on a 
S defined node set. 

00,1, ТУТ, Т 

КОЕТ NODE BEST 

1 

Wee 3, 47.3, 6, 1, 8 


625,626,627, 628 
5 Face Set wee 


61 


S 

*SET SEGMENT 

$ This is the DAA] boundary. 
2 

1236, 235 


927 6277 0255530 
* END 


2D USA input decks 


Of particular importance in 2D models utilizing the USA code is the numbering of 


nodes along the DAAI boundary. This boundary is located around the circumference of the 


fluid mesh. Imagining the 2D cross-section as a cylinder section, the first two nodes of 


these DAA1 face elements should progress around the circumference not along the X-Z 


axis. The NSORDR command in the USA FLUMAS input deck allows the user to 


cyclically reorder these nodes to achieve this requirement. This my have to be done 


visually and can be very time intensive. 


FLUMAS DATA FOR 2D BOX MODEL 


flunam geonam strnam daanam S FLUNAM GEONAM GRDNAM DAANAM 
T T T S PRIGMT PRTTRN PRTAMF CALCAM 
TAME Е S EIGMAF TWODIM HAFMOD QUAMOD 
F F T F $ PCHCDS NASTAM STOMAS STOINV 
FEET S FRWTFL FRWTGE FRWTGR FRESUR 
F T F F $ RENUMB STOGMT ROTGEO ROTOUA 
F F F F S PRTCOE STRMAS SPHERE ROTSYM 
F F F F $ OCTMOD CAVELU FRWTEV INTCAV 
F F F $ BOTREF MASREF BMWHIP 

0 628 0 52 S NSTRC NSTRF NGEN NGENF 
ООО S NBRA OAL NCAV 
0.9389E-4 59155.2 $ RHO CEE 

10 $ NVEC 

DNE UE TD. $ DEPTH CXFS CYFS CZFS 

14.7 386.088 S PATM GRAVAC 

0 $ NSRADI 

19 $ NSORDR 

DEUS ж! S NORD JBEG JEND JINC 

ЕІ S NORD JBEG JEND JINC 

оа S.NORD JBEG JEND ae 

LAS 1 $ NORD JBEG JEND JING 

028921 1 5 NORD JBEG JEND JINC 

ШЕ?) 22221 $ NORD JBEG JEND JINC 
002492571 $ NORD JBEG JEND JINC 

т S NORD JBEG JEND JINC 

028 29 1 $ NORD JBEG JEND JINC 

130957 9 SENORD JBEG JEND WING 

0 33 3597] $ NORD JBEG JEND JINC 





lc ЛА. > NOR Des 
Q D Sa $ NORD JB 
des Jae S NORD JB 
O 40 43 1 S NORD JB 
144 45 1 > NORD JB 
0.46 16M S NORD JB 
1 47 47 1 $ NORD JB 
а S NORD JE 
AUGMAT DATA FOR 2D BOX MODEL 

strnam flunam geonam prenam $ STRNAM 

Бие в Е S FRWTGE 

мал Вер ә ЕТПЕК 

EE EE 3 И РЕТСМШ 

ЖЕКЕ m $ MODTRN 

Il S NTYPDA 

628 1884 3 3 S NSTR 

1 5 INSBTEC 

В 52 1 рес 


IEMINT DATA FOR 2D BOX MODEL 


prenam posnam $ PRENAM POSNAM 
resnam $ RESNAM WRTNAM 
EE E $ REFSEC FLUMEM PWACAV ITERAT 
DE Э ӘСІШСЕТЕ СЕМТЫП BUOYAN 
10 $ МТТМТ NCHGAL 
ШІ) 1.05-5 а ЛЕГ DELTIM 
мВ ЕЕ $ EXPWAV SPLINE VARLIN PACKET 
ШИШЕ Е SOHYPERBOEXPDOS DOUPDC S Fp pi 
EE Е Е > BUBPUL ЭНЕЕОВ 
1 5 МСНАКС 
on зи DEPRE 
E 0, -132.0,6.0 SES YG Le 
mee. O,-40.0,6.0 S от SZ 
201 SUE ESTE 
O. $ PNORM DETIM 
6.4E-6 S5 DIHIST 
2 2 СПСЕСР 
En 206.0 25.0 $ WEIGHT SLANT CHGDEP 
2000 2000 SONSAVERSNEESET 
оо О S LOCBEG LOCRES LOCWRT N TARI 
Е Е F F S FORWRT STBDA2 ASCWRT 
EL .0,-16.0,6.0 S XV YV ZV 
E SHETSPLZ 
3D LS-DYNA input deck 
* KEYWORD 
ЕТЕ 


BARGE W/OUT BEAMS 
*CONTROL_ TERMINATION 
Pree. ,0,070,.0 
“CONTROL TIMESTEP 


EG JEND 
EG JEND 
EG JEND 
EG JEND 
EG JEND 
EG JEND 
EG JEND 
EG JEND 


ЕТОМАМ 
FRWTST 
DAAFRM 
РЕ ТЕН 
TRECE 


NSFR 


JSTART 


JUNE 
JINC 
JINC 
ЕС 
ЈІМС 
ЈІМС 
ЈІМС 
JONC 


GEONAM 
FRWTFL 
SYMCON 
PRESTE 
INTWAT 


NFRE 


JOP 


PRENAM 
LUMPFM 
DOFTAB 
PRTAUG 
CEAPRE 


NE TR 


TINC 


Е 000060 
*CONTROL PARALLEL 

Zo ul 

*DEFINE CURVE 

| 

0.,1.0E-5 

0505, 1.0533 

*DATABASE HISTORY NODE 
ies, 1 65 596,005 96 
AO 2853, 267 

* DATABASE NODOUT 

000002 

*DATABASE HISTORY SOLID 
2257, А200 БОТУ МЕР ат 
22500 (413253222, 6327 

9557 95329706342 

*DATABASE ELOUT 

0.00001 

*DATABASE BINARY D3PLOT 
0.001 

*DATABASE BINARY D3THDT 
02003 

*DATABASE EXTENT BINARY 
(UMS Logs l.1,1 
OP. 70. 0. 0 

*BOUNDARY USA SURFACE 
а-ы) 

РТ Ь DETONATION 

Е О, =132.0.-2898. 000 
1, 0. 000321,84.0,-16.699,-4.0,. 3721 
° 

*NODE 
7202.000000000E+00,-12.0000000,02000000000E 007970 
2321184. 0000000, =-652.0000000,-300.000000 7070 
*PART 

púa (acoustic) 

90,74, 90 

SECTION SOLTD 

4,8 

*ELEMENT SOLID 
aI ALG, AA, EO 


Ф о е о ә 


20592790, 22439,22171;22174,;,22491,23266,;, 2290 123268 
“PART 

Shell Plating 

Soo E 

*SECTION SHELL 

Se pees 
ООО 25, 0225 
ELEMENT SHELL 
о 

$ 

В 
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5 
S Fa 
= 
SE: 
1 


ce set wet 


SEGMENT 


699,698, 697/696,0.000Е-00, О  ИОК-ФОФ РОзОЕГООО ОООКТИВ 


771557720 ,8497,8492, 0.0005 00 0500 OE 00] ЕО, 0. 0005306 


*MAT ELASTIC 


1,7.350E-04,3.000E407, 0.300 


5 


I T ACOUSTIC 
09455-05, 5-9Т6Е 0470-5018 засы 


ШЕШЕГІ 070 0. 1.00 


*ELEMENT_MASS 


1,69 
El 
СЗЭ 


‚0.664 
5,0 2604 
0,0.664 


* END 


FLUMAS DATA FOR 3D BOX MODEL 
flunam geonam strnam daanam 


E 


S 


нат fm] 


3 


21 
4.7 


C O t2 pF UJ 0 C) СО ir VW rrj 


AUGMAT DATA FOR 3D BOX MODEL 
strnam flunam geonam prenam 


в Е 


pon nj nj 
"nj 3 
n] nj nj 


1 
2,21 
1 

Ue 


3D USA input decks 


Е 


nj nj nj ej nj nj d 


ES 
T 
Е 
Е 
Е 
Е 
Е 
ОТТО СЕ 

0 

45.249591 008 


x. 
386.088 


bb 
Т 
а 
Е 


7698810333 


3906501 


AP ФФ Ехо 


+» Ар Ар Ак +) Xp n Xp «Xn» 4p -U X0 Xp 40 X0 X0» 


FLUNAM 
ERTEMT 
EIGMAF 
PCIICDS 
FRWTFEL 
RENUMB 
ERTESE 
OCTMOD 
BOTRER 
NSTRE 
NBRA 
RHO 
NVEC 


GEONAM 
PRTTRN 
TWODIM 
NASTAM 
FRWTGE 
STOGMT 
STRMAS 
CAVFLU 
MASREF 
NSTRE 
NCTE 
CEE 


GRDNAM 
ЕКТАШЕ 
HAFMOD 
STOMAS 
FRWTGR 
ROTGEO 
SPHERE 
ЕКИТЕУ 


NGEN 
NCAV 


DESTE CLES CIES CEES 
PATM GRAVAC 


NSRADI 
NSORDR 


STRNAM 


FLUNAM 


GEONAM 


DAANAM 
CALCAM 
QUAMOD 
STOINV 
FRESUR 
ROTQUA 
ROTSYM 
INTCAV 


NGENF 


PRENAM 


FRWTGE 
FDUSEY 
PRTGMT 
MODTRN 
NTYPDA 
NSTR 

NSETLC 
NDICOS 


ERES T 
DAAFRM 
ERTEERN 
SIRECL 


NS ER 


JSTART 


ЕКИТЕТ 
SYMCON 
PRISTE 
INTWAT 


NERE 


JSTOR 


LUMPFM 
DOFTAB 
PRTAUG 
СЕЛБЕЕ 


NFTR 


JINC 


TIMINT DATA FOR 3D BOX MODEL 


prenam posnam 
resnam 
ER E 


tr] 


îj 3 Hr] O O f] 


Q p Hj rj Hj O m Hj "1 
mj =} j] = 


SE O NL SL c2598.0 
83.0, 64:0, =-00.0 


201 

Too 

6.4E-6 

2 

СО 07250 
2000 2000 
0000 

EE Г 


Р.О 


la 


O NORNER N ARN ЧУУ Y AENA AH AN ANAND AN AN 4л N xu 


PRENAM POSNAM 
RESNAM WRTNAM 
ПЕББЕС ELUMEM 
INCSTRASENINT 
NTINT NCHGAL 
ЪЛЕТ ТИ» ВЕСТ ТМ 
EXPWAV SPLINE 
П СВЕЦЕ Е ДЕШЕ 5 
BUBPUL SHKBUB 
NCHARG 

HrDERE 

XC LG 

SX SY 
ШЕБІ 
PNORM 
DT! T5SE 
CHOTYPE 


DETIM 


PWACAV ITERAT 
BUOYAN 


VARLIN PACKET 
DOUBDE VELBLINP 


Ze 
SZ 


WEIGHT SLANT CHGDEP 


NSAVER NRESET 
LOCBEG LOCRES 


LOCWRT NSTART 


FORWRT STBDA2 ASCWRT 


XV IV 
DISPLA 
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ZV 


APPENDIX C. FLUID MODELING USING TRUEGRID 


This appendix covers the procedure for creating a fluid finite element mesh using 
TrueGrid. Some of these techniques were adopted from Reference 17. Two of the most 
useful features of TrueGrid for coupled fluid modeling are the BLUDE and BB commands. 
The basics of using TrueGrid will not be covered here and some familiarity with the code is 
assumed. An example of a partial batch file is Included. This batch file creates a basic 
barge and extrudes a fluid mesh. Additional information can be found in the TrueGrid user 
manual. (Ref. 8] 

The BLUDE command "extrudes" the structural mesh through a user defined grid 
or block of mesh mated to the structural wetted surface. The block part 1s actually attached 
to a surface definition created from a faceset of the wetted elements of the structural mesh. 
Тһе resulting extruded mesh matches exactly to the structural mesh, a prerequisite for 
successful fluid modeling. 

The extrusion procedure is as follows, with important commands and menu 
selections denoted in bold and all capital letters for emphasis: 

1. A structural model must be created. TrueGrid can be used or the READMESH 
command can be used to input a mesh from another code format. It is very 
important to remember though, that when TrueGrid reads in a finite element mesh 
from an outside code format, it renumbers every element and grid point. Therefore, 
once the mesh is through being manipulated in TrueGrid, and it is written an output 
file, the grid point and element 1D numbers will not match between the original and 
newly output model from TrueGrid. 

2. The elements of the structural model that will be in contact with the fluid, 1.e. 

the wetted surface, must be grouped into FACESETS. This option can be accessed 

from the environment window under the PICK option by choosing the SETS 
button. The FACES button should be selected. Appropriate faces of the models wet 
surface should be put in a separate FACESET. Pick faces which are naturally 
defined by the geometry of the wet surface. For a rectangular barge this includes 
the bottom, sides, bow and stern below the waterline. For a ships hull this would 
include the port and starboard sides and the stern. In this case the bow is typically a 
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sharp edge and would not be selected as a FACESET. The HIDE drawing mode 
vice WIREFRAME should be used for the mesh to ensure that only the visible 
elements are picked. This will make faceset selection must easier, since it must be 
done by hand using the lasso tool guided by the mouse. The three- or four-node 
selection option is the best to use when choosing the faceset. This means that three 
or four nodes of an element must be within the selection lasso for the element to be 
added to the faceset. The selected elements will be highlight in white. If some 
elements are selected that are not desired in the particular set, they can be easily 
selected and removed; using the one node selection option is best for this operation. 
The REMOVE button should be pushed also. The set must be named and saved 
once selected. 

3. The SURFACE menu SD (surface definition) option should be chosen next. A 
surface number must be input. The faceset option should be selected from the end of 
the surface options list and the name of the desired faceset should then be input. 
This step converts the named faceset into a surface definition. The new surface will 
be displayed in red in the physical window. This is also useful in troubleshooting 
the faceset. The SURFACE created should have no holes in it. These holes would 
represent elements that were missed in the FACESET selection. 

4. Next, the PARTS menu should be selected and the BLUDE option chosen. 
Using this option, the user creates a block part that will be attached to the surface 
created above. This block will serve as the "guide" for the extrusion of the structural 
mesh; therefore, the block's mesh must match the structural mesh or be of finer 
quality in order to get a quality extrusion; an exact match is not required however. 
This block part is created in the same way as a block using the BLOCK command. 
The blude command requires two additional inputs, however. First, the name of the 
face of the block where the extrusion begins must be input. This is simply the face 
closest to the structure. Next, the name of the faceset to be extruded must input. 

5. The block part created can now be manipulated to obtain the desired geometry. 
This is the “art” behind TrueGrid and requires practice to be done successfully. The 


part created must now be attached to the surface defined in step 3. It can be attached 
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using any of TrueGrid's available options. The PROJECT button in the 
environment window 1s the easiest to use. This will work for simple cases, but a 
complex surface may require use of other TrueGrid methods. 

6. The interface of the extrusion mesh and the structural mesh should be carefully 
examined. Orthogonality of the fluid and structural mesh is a must (next to the 
wetted surface) and should be verified; TrueGrids DIAGNOSTICS menu provides 
the necessary tools. The block mesh can be modified as needed using various 
TrueGrid tools to ensure a quality mesh is constructed for the extrusion; two 
examples of useful tools are the mesh relaxation algorithms and use of a cubic spline 
to added curvature to the block mesh edges. Material properties can be assigned to 
the mesh also, just as with any other part in TrueGrid. 

7. Once the user is satisfied with the mesh, the MERGE command should be used 
to end the PARTS phase. The result will be a fluid mesh, which matches exactly to 
the structural mesh. The mesh will consist of 8-noded solid elements. The STP 
option should be reviewed and used to fully merge all the nodes on the interface of 
the structure and the extruded mesh. Prior to merging, the extrusion mesh can be 
replicated using the LCT and LREP commands. This will only be effective if the 
model is symmetric. Using these part replication features, the user only has to build 
one-half of the extrusion mesh. 

8. Additional extrusions can be performed, including on any newly extruded mesh 
surfaces. This must usually be done to fully form a fluid mesh around the structural 


model. 


The BB command defines a block boundary interface for parts that will interact with 


one another. This command allows the user to save the geometry of a face, edge or vertex 


of a part and use it to form the geometry of another at a later time. This feature 1s used in 
the PART (or BLUDE) phase. The useful features of this command for fluid modeling will 
be highlighted here. 


l. In the PART phase the user selects the face of a part for future reference and 


gives it a specific name. This is done with the BB command, the region and a name. The 
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user can select the region in the computational window and print the region with the Е! 
button. For example BB 11155 1 12, defines the master side of a block boundary as 
region 1 1 1 25 5 1 and names it 12. 

2. Use ENDPART or MERGE to complete the above part. 

3. In the PART phase the user selects the face of the part that must match the 
MASTER defined above. Use the same technique. Now define the BB command again 
with the newly selected region and the same block boundary number as the MASTER. For 
example BB 112222 12, matches the region 1 12222 to theregion1 115 5 I. 


Example of a TrueGrid batch file to create a coupled barge/fluid model. 


Material Definitions use the mate command. 


c mate 3 is structural steel for shells 

ETNIES 3 l1 sbruct rho ./35e-3 e .3009209 pr 2300 + 

€ mate ОО s tluid 

и SGU SO THO 22836-44 SS 59155 2 D 59 ct 1 Езр Т ІІ сет Ale 


The block command creates the basic barge structure. The negative signs generate shells. 
The thic command defines the shell thickness. The mate command assigns a part and 
material number. This will be the part number used in the LS-DYNA input deck. 


block -1 -6 -12 -17; -1 4 -7; -1 7; 
0750 118 EORNM =12 0127 0 24; 
ENTE. 25S Male 


merge 


Use the pick sets command to pick the bottom faces of the barge. Defıne this set as a 
surface with the sd command. Do the same for the bow and stern below the waterline. 


fset bottom = ls 

c linear shells 

1963232096213 1198221722347 31923 5,0 352.065: 
sd 1 faceset bottom 


[Sc EMEN S 
c linear shells 
1 :: 


fset aft = ls 

c linear shells 
3924056415. а: 

sd 2 faceset fwd 
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sd 3 faceset aft 


Use the blude command to generate a block of mesh to extrude the faceset on the bow of 
the barge through. Do the same for the stern. Define the fluid mesh as mate 90. The res 
command is used to vary the mesh density in the fluid volume. 


blude 2 fwd 

1411; 1.278122, BALO o ua 
res 1 1472 2 27939995957 

mate 90 

u 2,0208 

merge 


Buades” 1 aft. 1 11; 1.7;.1 4;168 2522 -=12 120 12. mate 90 
mes) Las ee, 

сле Ба 3 

merge 


The above commands generate two sections of the coupled fluid volume off the bow and 
stern. Next two facesets are defined along the port and starboard sides of the new coupled 
barge/fluid model. 


fset stbd = 153 
c linear bricks - face #3 
EM 6078190 ls Ome dee 50a 8 1 2905341: 350. 
c linear shells 
И. 


sd 4 faceset stbd 
eset port = 155 
с linear bricks ~ face #5 
щат ЗОРО: ope e le yO 201:210 53521: 330; 
c linear shells 
I 0-55 145-1752 289303: 


sd 5 faceset port 
The fluid is extruded out to the port and starboard with the next commands. 


Exude 4 port 1 60;1 8;1 4; 
Оо 252;-64 -12;0 12; 

mate 90 

Nes 11 1 2=2 2 1 -96957 
eee, —2;;sd 5 


merge 


meade S ~stbd 1 60; 1.67 1) 4: 
K 1 252;12 64:00 12; 

mate 90 

O do2 2-2 15 

SE ; -І;;за 4 

merge 
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The next command defines a faceset on the bottom of the coupled barge/fluid model. This 
is then sxtruded to the desired depth. 


fset bottom = 1b2 
с linear Dri cks -~ face #2 
о 02D UOS SIId0SIS]-700 211:220-281:290 3115320 3515360 


сте АО УЛС 775: 1791 13963 3 09 1024:1080 1088:1814 18 
c linear snells 
ЕОС ОЭ оаа LOGE SLI SS CS LOOL 
sa 6 faceset bottom 


blude 6 bottom PESO 1 21; 9042925205864 54; —BN2090 90 


J 
ра 
sí1 722225076 
merge 
RP 1 


Now that the entire fluid volume has be created, a faceset is defined on all of its exposed 
faces except the free surface. This faceset is the DAA1 boundary. The orpt command 
ensures the nodal numbering of this faceset set generates outward facing normals. This is 
crucial for proper behavior of the analysis. A good trouble shooting technique is to create a 
surface of this faceset (sd 7 faceset wet) and then remove the mesh with the rap command. 
This surface should no holes in it and on the edge along the free surface should be visible 
when it is rotated. The dap command is used to restore the mesh. 


eub - 0.00 

fset wet = lbl 

C linear bricks - face #1 
ООО -0 60 70 80 90 


sd 7 faceset wet 


The output command can now be used to generate an LS-DYNA input deck. 


APPENDIX D: SECOND GENERATION DDG53 COUPLED MODEL 


This Appendix details the modeling of a second-generation coupled structural-fluid 
model for the simulation of the shock trial of USS JOHN PAUL JONES (DDG53). Some 
of the results of a simulation conducted using a first-generation model were presented in 
Ref. 7. This paper showed promising results in simulating the vertical response of the 
DDGS53 when compared to the ship shock trial data. Further investigation into the coupled 
structure-fluid model was desired to improve the nodal distribution in the coupled fluid 
mesh. Specific areas of interest were along the keel line and bow. The extreme geometries 
of the sonar dome presented a significant problem in the meshing of the model. This model 
uses the same three-dimensional finite-element model of the USS JOHN PAUL JONES 
(DDG53) developed by Gibbs and Cox, Inc. [Ref. 22]. Table Dl compares the 


characteristics of the old and new coupled models. 


Table D1: Comparison of Coupled Fluid Models 
Characteristic Old Model New Model 
Number of Total Nodes 174,635 198,130 
Number of Fluid Nodes 139,664 162,930 
Number of Ship Nodes 34,971 35,200 


Number of Fluid Elements 94,199 132,810 
Number of DAA] Faces 4,290 5202 
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Structural hull sm Í 
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Figure D1: Sketch of Modeling Scheme 


Figure D1 shows the basic scheme utilized in constructing the fluid volume. This 
scheme seems very effective in dealing with the complex geometries of a ships hull. It 
shows the structural hull fitted within a single element thick liner. This liner is used to 
smooth the discontinuities in the ships wet surface. An inner liner is then extruded from 
this single element liner to give the fluid volume a cylindrical shape. The inner cylinder is 
fitted into an outer cylinder that represents the bulk of the fluid volume. Bow and stern 
caps are then extruded from the inner cylinder to plug the ends of the outer cylinder. 

The intent of this Appendix is to provide a near step-by-step approach to 
constructing a coupled fluid volume to a ships hull. The DDG53 model will be used as an 
example of the procedure. The courier 10 text is from the actual TrueGrid batch file. It is 
included here as an example of the command syntax. This is NOT a complete batch input 
deck. This model was constructed with the assistance of Dr. Robert Rainsberger of XYZ 


Scientific, Inc. His help was invaluable. 


c this file loads the ddg nastran deck and deletes the springs 


parameters or 1344 c outer radius 
bh 1000 c distance behind the stern 
ООО ШШ ce smoothing flag 
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The readmesh command reads in a supported codes mesh and converts it to a TrueGrid 
part. The endpart command is used to finish that part definition. 


readmesh nastran sp 2nfÜf.bdf 
endpart 


merge 
The delspds command deletes the defined numbered springs. 
delsods 15421 6; 


The block command 1s used here to create a patch of mesh to plug the hole in the ships hull 
where the propeller shafts extend through. 


€ fill the shaft hole in the ship 
БШӘӨСЕ 1 2234 56 7 8; 


=1; 

I a 

1152 120041228 1295 1344 13947 1440 1488 
525 

o0 7l 8288 


endpart 


c fill the shaft hole in the ship 
mock 1-2 3 4 996 7; 
-1; 
ШЕТА; 
1152 1200 71:249527200 884045] 3920] 40 
135 
GO 7132279 


endpart 


merge 


The pick sets option, described in Appendix C, is used here to create two faceset’s of the 
wet surface of the ships hull, starbrd and port. This translates to fset in the batch file. Two 
surfaces are then created from these facesets, sd 1 and sd 2. 


C create a face set of all of the polygons forming the exterior of the 
ship 


fset starbrd = ls 
fset port = ls 


C form the two surfaces of the wet ship 
за 1 faceset starbrd 
за 2 faceset port 


C extract the edges of the new surfaces as curves 


cura | зе 1.2 
cura Z se 1.3 


тӘ 


cura 3 se 2.1 
curd 4 se 2.3 
с create curves that form the special region along the keel 


The pick sets option is used here to create a faceset of the bottom of the sonar dome. 


c form a face set of the bottom polygons along the sonar dome 
o 15 ]223:1260;; 


The blude, described in Appendix C, command is used here extrude the sonar dome mesh 
one element orthogonally. 


С 
blude 6 sonar I 61:1 20 39:7 2;4896 55812 58 -9 0; 


endpart 


The block and bb commands are used here to generate a block boundary interface. The bb 
interface will be used to form the single element liner. 


c create the starboard block boundary interface part 
c the part is not used, just the block boundary interfaces 
ОВК 19а О 


DSL 
5S 1 2 
So ie 


Nie ва 
N N N 
N N аз 
Юон № 
е 


` 


mate 0 


endpart 
The blude command is used here to extrude the starboard side mesh one element 


orthogonally. 


с create the first orthogonal layer of elements along the starbooard 
ае starnord 1 194.1 2.1 24;0 5584 99050 254 


Ори dz ¿1 : 

реа ie T Normals: 
peel pina P s s 517.153, 

ЕЕ A 57 

endpart 


The block and bb commands are used here to generate a block boundary interface. The bb 
interface will be used to form the single element liner. 


c create the port block boundary interface part 
c the part is not used, just the block boundary interfaces 
Sie qo c —IX21.0 95900 0-0 7254 


bb 2 жо 
bb 2 1 4; 


N P 


s wa 
TCU 
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bb l2 SHE 

пате 0 

endpart 
The blude command is used here to extrude the port side mesh one element orthogonally. 


c create the first orthogonal layer of elements along the port 
blude 3 port 1.194;1 2;1 24:0 5500. 000 51 


bb 
bb 
bb 
bb 


; 
normal -9; 


, 


eere нњ 
UNNE 
N H| ҥн P> 
NNN N 
N N Ñ P 
Г в КО КО 
бу P КО КО 


endpart 
merge 


The port and starboard liners created above leave a gap along the keel line. The block 
command is used to create a series of wedges to fill in this gap. 


C Create the wedge elements filling in the gap 
block 1 2 3 4 5 6 7 8 2 
TOTAL Жы 22-222 2502 12 70015573425 
20.21 228226 AS oo 28825 
STD ITT 33627 38 39 
40 41 42 43 44 45 46 47 48 49 
503651 52% 59: 50 ФЕ 57-55. 759 
СОМЕ 22055020277 55 22 GT 68 69 
0-71 2 ста И ES, 
ӘП” 8102782 бес 29200200050 000 
Do Е O ME 
100 5012102810 91019 05 106 107 1931103 
POE ТР 15 Ole zer 
ШОШ ho] c 123712 2171125 Об ТОТО О 
1302151. 1232. 1357; 
ор 


| 


SO DIS DOTES 
о 
о 
O O O OOO 
оо 
aa со CO €) 
OOO OO 
СОЕ СЭСО С O 
QC O CO © 
соосу 
(DSS SS, 
SOS] DAS 
OO OOS 
SOS ав 
OOS © 
e o a 
о о Со 
aaa a 
EI) 
OSS SOS 
OSO OO © 
SS OOS 
SOS © 
C» SOE 
ODE 
SOPAS C) 

(® ®) => @; 
OOO 
OOO O 
COITO 


The dei command is used to delete elements. 


II E A 


endpart 
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The sd and cy commands are used to create two surfaces that will define the inner and 
outer fluid mesh cylinders (see figure D1). 


Sd ll Gy 0 0 254 ] 0 0 529 c intermediate cylinder 
замест с OD AI GO USO cC cuter cylinder 


The pick sets and sd options are used to create a surface of the stern. 
fset stern = 1b2 


sd 23 faceset stern 
tee 5554 521720354569 7 80971002112 


endpart 
The block command is used to create a block of mesh on the stern. 


pillo ci 2102 34 -5/182.>4;5.5513301e+09 455501850] e+ 038-070. 0 CASOS 
се 1 2 0 4 Sal 2 70074; 


merge 


The pick sets and sd options are used to create a surface of the exterior faces of the single 
element liner. This faceset and surface will be used to extrude the bulk of the coupled fluid 
volume. 


c make a face set of the exterior faces of the first orthogonal layer 
fset pull = 101 


sd 10 faceset null 
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Original Mesh 
Mesh after Initial Manipulation Ruid Mesh after projection onto 
Prior to projection and Extrusion Hull, inner and outer cylinders 


Figure D2: Extrusion of Bulk Fluid 


The blude command is used to extrude the bulk of the fluid from the faceset defined above, 
fset hull. This block of mesh 1s manipulated to form the desired geometry for the extrusion. 
The mesh is then projected onto the inner and out cylinders defined above. The bb 
command is used here to create a set of master block boundary interfaces to be used for the 
bow and stern caps. See figure D2. 


Elude 6 NULL 2 0 oclo 
ото. 
>; 
=1.3062277e+00 8.6497479e+02 9.624691772+02 17056408 Se 05 
4.2724814e+03 4.9031758e+03 5.1360000e+03 5.4480044e+03 
1 5921299е+03; 
-300 -90 -87 0 87 90 300; 
= (O) -408 0; 
msec i 28 2 2 80 22 7 10 8 
mseq 195035950019 
mseq ВТБ 


Ser 


Do 9 
сее 
po l 
BP 1 


EXPE 
= = [=> = 
= = og 
м ل‎ Jos 
N N N N 
F2 4o 0 -J 
OE 


endpart 


The pick sets command is used to select the fantail of the ship model. This set is then 
extruded to form a block of mesh. 


жые пап = 15 
Poo aces et pan 


19 


в гр. Тап 1 lla 17;1 6; 
urn 20,223. 228 75;219.5 304; 
ee oC) | 


The pick sets command is used to select the faces of the stern fluid needed to construct the 
stern cap. The blude command is used to extrude the mesh of the inner cylinder to the end 
of the model. The bb command is used to define the slave side of the block boundary 
interfaces, 9 and 10, defined above. This allows the mesh that is extruded from the inner 
cylinder to be projected down onto the outer cylinder mesh. The btol command is used to 
merge the nodes between parts 12 and 14. Part 12 1s the bulk fluid and part 14 is the stern 
cap. A similar procedure is used for the bow cap. 


1102 ا ا )ته 
В < prt strn‏ 
Eset Stb Stringer. 1b3‏ 
о 5п‏ 
fset SM, flu = 1b2‏ 
Semeemetaecset Strn flu‏ 


Ede 2 strn flu 1 lll l p S. 8 В. 75 0 228.75; 
-156 219; 


ВВ 1221 10; 
a 1 223 179; 
merge 

ВЕЩИ 


БОЕО Т5 144 
БЕСОП 12 ld 


The procedure for the bow cap 1s similar to the stern cap above. 


с=с plug u= пра 
sd 42 faceset fwd plug 


БЕС рот = 162 
sd 43 faceset fwd port 


Во - 
Ве 
DIGE RETA BIG LIL 22 45.1 21, -904 60:02, 


. 
, 


bb 7 
b 8; 


Tee 
О Е. 
bptol 12 15 3 


The fluid mesh is complete. As detailed in Appendix C, the outward faces of the completed 
fluid mesh are selected with pick faces command. This faceset 15 the DAA1 boundary for 
the model. 
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The following figures show the details of the modeling of the USS John Paul Jones 
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Single Element L 


Figure D4 
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Complete Single Element Liner 
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Figure D6 
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Figure D9: Coupled Ship and Fluid 
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Figure D10: Top View of Coupled Model 
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APPENDIX E: DAA1 ON WET SURFACE/NO FLUID VOLUME 


For completeness a comparison is provided of the response of the barge model with 
and without the coupled fluid volume. Figure El shows the vertical velocity response of 
node 159 for the coupled and uncoupled cases. As can be seen in figure E2, the initial 
response at this node to the incident shock wave is similar. Both rise rapidly to 
approximately 25 ft/sec although the uncoupled model lags the coupled model slightly. The 
response curves then separate with the uncoupled model quickly dropping off to a constant 
positive velocity of approximately 2 ft/sec at 5-msec. The effects of cavitation are seen in 
the coupled model response as a higher peak velocity, longer decay time and larger 


amplitude rider frequency. 


Coupled vs. Uncoupled Response for Node 159 


— Coupled Model ——No Fluid/Uncoupled 





Vertical Velocity (ft/sec) 





Time (msec) 


Figure E1: coupled vs. Uncoupled Response for Node 159 
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Coupled vs. Uncoupled Response for Node 159 
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Figure E2: Initial 4-msec of Response of Node 159 


Figure E3 shows the comparison of the models at node 213. This 1s the midships, 
centerline node. The response is damped at this location due to the flexibility of the plate 
and the presence of a lumped nodal mass. 


Coupled vs. Uncoupled Response for Node 213 
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Figure E3: Coupled vs. Uncoupled Response of Node 213 
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